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РНЕКАСЕ 


- 

Precipitation static has long been recognized as a particularly serious 
problem in aviation radio navigation and communication systems. Many investi- 
gations as to its causes and nature have been иас u with successful 
progress, especially during the past fifteen years. Static free navigational 
systems (VOR Omni-range), for short ranges, and line of sight communications 
systems (VHF-UHF), are now standard equipment on all aircraft, although it 
must be pointed out that these solutions involve a sort of strategic re- 
treat, that is, a shift to higher frequency levels where the precipitation 
static energy content is of a lesser value. Also the investigations that 
were made all utilized propeller driven aircraft for their test vehicles 
and although projects were earmarked to include jet type of aircraft, none, 
to the author's knowledge, were ever reported on. 

The problem that is the subject of this thesis was essentially to de- 
sign a device to measure the relative noise power or voltage due to corona 
discharge, at certain specified locations in the frequency spectrum, and 
be carried aboard jet type aircraft. This is a bit of a swift expianation 
but a detailed account will be given later. 

Also an attempt is made to correlate information regarding some of 
the statistical aspects of noise analysis and altnough it is not a completely 
rigid treatment, it describes in a condensed manner, a description of noise 
signals and their effects on electronics circuits. 

The author would like to express great appreciation to Messrs. Н. Р. 


Blanchard and Maurice Mills of the Stanford Research Institute for 
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Chapter I 


INTRODUCTION 


l. History and General Considerations of the Problem 

Precipitation static may be described as radio interference caused by 
corona ; It is sometimes labled as, 'Poisson noise', is tech- 
nically named, ‘impulse noise!', and is statistically described as, 'nearly 
normal random noise'. It is but a single type of the many kinds of noise 
that have been studied, analyzed and catalogued, although it does have the 
dubious distinction of being the number one culprit of aviation naviga- 
tional and communications systems. | 

The first analysis and explanation of precipitation static was that 
charged particles, created by the turbulence and charged energy of thunder 
clouds, generated static when striking the antenna. No attempts were made 
to sie ds these effects for many years until the introduotion of the 
shielded loop antenna in 193h. It was assumed that since particles did not 
strike the antenna proper, static would be reduced. However, this assump- 
tion was not substantiated in practice because it was found that in actual 
air operation, severe static easily disabled the system. 

In 1939 after more investigation, a new theory 11 assigned the main 


cause of precipitation static in aircraft to that of corona discharge. 


Later the concentrated effort of the Army-Navy Precipitation Static Project 8 


verified this theory which resulted in the subsequent changes to air-frame 
construction and airborne electronic equipment that are in standard prac- 
tice today. 


with the increased use of jet and rocket propulsion in supersonic 
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aircraft and guided missiles, it was erroneously assumed that the determi- 

nental effects due to precipitation static would be of no consequence, 

since in this type of craft, a natural discharging process due to the ther- 

mal ionization of the exhaust gases, should efficiently transfer free charge 

to the atmosphere. This 18 partially accomplished because a complete recom- 

bination of the ions produced in burning the fuel has not as yet taken place, 

thus the electric field of the aircraft will capture those ions of opposite 

Sign while simultaneously repelling those of the same а 
In a recent report by R. Г. Таппег,201% маз зћомп that the strength 

of the noise signals due to corona discharge is directly related to the 

dipole moment of the discharge current. The dipole moment is proportional 

to the product of the current flow and the distance over which it takes 

place. Thus in any supersonic aircraft which strikes ice crystals (cumulo- 

cirrus type clouds), or heavy turbulence (cumulo-nimbus type clouds), or 

the layers of ionized air located at extremely high altitudes, a static charge 

builds up in the vicinity of plastic antenna housings, canopies and control 

surfaces. It is here that one is faced with a new type phenomenon that is 

based on a very old principle: that substances of high dielectric constant 

acquire a positive charge when rubbed with substances of lower dielectric 

constant. This is called triboelectric charging, and can be visualized as 

a scrapping off of the electrons by the airstream. Naturally this charging 

eventually culminates in a corona type discharge at many points close to 

the outer perimeter of plastic housings and the dipole moment of the cur- 

rent discharge acts as an efficient antenna to couple this electrical 

noise into the receiver input circuits. 


The problem specifically, is to determine the average noise power or 
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tic, recorded for a period of 
time equal to the flight duration of the aircraft. This system must dupli- 
cate an actual navigational system's characteristics with respect to band- 
width, resonant frequency, servo system time constants and antenna match- 
ing stage. 

The actual navigational и сап be represented in general by the 


diagram in Figure 1. 
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FIGURE 1 


Simplified Diagram of Navigational System 








To illustrate: 

If the incoming waveform is sinusoidal in nature, then the output of 
the mixture will be the sum and difference terms of the original inputs. 
This is applied to the integrator which acts to average only the difference 
terms since it is a low pass filter. The output will then be a very slowly 
varying voltage and the subsequent error signals will also be very small 
in magnitude. 

Where precipitation static is present the waveforms (idealized) might 
take the form of those shown in Figure 2. Since power is proportional to 
the voltage applied to the servo system, the disturbance to the system is 
related to all the input noise pulses integrated over a period of time 


equal to the time constant of the system, and can be written as 


«x 

\ 

T E dt = Phase Disturbance in System 
о 


The error signal to the reinsertion oscillator will then be proportional 
to the input average noise power so that periods of heavy static would pre- 
cede large error signals. 

The system that was designed and constructed duplicating the above 
characteristics is now being operationally used to record the average noise 
power over periods of ten to twelve hours. These results are later com- 
pared to the navigational errors plotted separately by the navigator to see 
what relative values of average noise power cause the maximum errors 18 


the system. 
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A second system to measure the average noise power at specific fre- 


quencies is alse described in Chapter III with some data which shows re- 


sults of actual flights made using the systen. 








2. The Projected Method of Solution: 

A study of the available systems for measuring noise due to atmospheric 
sources and corona discharge indicates a common tendency to make the major- 
ity of measurements using an electric field intensity meter, or circuits of 
some such similar nature. The results are then interpreted using a graph- 
ical portrayal of corona discharge current as a function of electrical field 
intensity. Іп this problem it is desired to measure the energy components 
of the noise at a specified frequency and known bandwidth, therefore a 
relationship with respect to time must be set up to measure this noise power. 
But the instantaneous noise voltage induced in the antenna cannot be con- 
sidered as a regular function of time, for if it were, both amplitude and 
phase relationships would have to be known for one complete period. Since 
it cannot, the problem is one that is statistical in nature. 

Statistics may be succinctly defined as the science of reduction and 
analysis of data. Its application to the problems associated with com- 
munication processes can be mainly attributed to N. Wiener$? who was one 
of the foremost contributors during the development years of statistical 
prediction theory, particularly concerning stationary random time func- 
tions, a random time function being defined as a function which does not 
follow any specific law of variation and tnerefore is not precisely predic- 
table. The type of noise being investigated in this paper falls into this 
category for it consists of a superposition of impulses occurring at ran- 
dom times with overlapping individual pulses. It is sometime referred to 
as Poisson noise because the Poisson distribution (Appendix I) is used to 
describe its statistical structure. In statistical applications the 


Poisson Distribution often appears when concemed with the number of 





occurrances of a certain event for a very large number of observations, 
where probability for the event to occur in each observation 1s very small.” 
It is for this reason that it is sometimes called the "small probability law". 
For the case of very high density Poisson noise, that Үл, when the number 
of elementary random impulses is large, and their mean duration is not 

too brief, a heavy overlapping occurs and the result is nearly normal ran- 
dom noise. This, according to Middleton?! can be considered to have a 
normal distribution with negligible correction terms, and the shape of the 
pulse is unimportant. (Appendix I) | 

It is plausible to assume that the corona discharge that takes place 
will have the characteristics of high density Poisson noise, for when the 
same navigational system was operated in slower propeller driven types of 
aircraft, the system suffered erratic operation (intolerable errors), 
only when thunder storms set up exceedingly heavy corona conditions. Also, 
the system operated normally under conditions of moderate corona indicating 
that the low density types can be ignored. 

For the case of the jet aircraft, the corona discharge intensity would 
be mainly a function of speed and since the latter is usually close to 
maximum designed speed, the intensity will always be close to maximun, 
with corona discharge pulses occurring at all of the places where plastic 
surfaces are exposed to the airstream. The mathematical theory to des- 
cribe this phenomena is covered below and is closely allied to the more 


familiar treatment of fluctuation noise. 


3. А development of Mathematical Description of Random Noise. 

Consider a plane-electrode diode for the emission limited case. A 
graph showing the current versus time for a single electron leaving the 
point of zero potential ( t = 0 ), and arriving at the point of highest 
potential ( t = t,), is shown in Figure 3. When a constant potential is 
assumed over a time interval 2T, and a pulse having the shape of that 
shown in Figure 3, takes place for each electron crossing, a plot of 
the current I(t) versus time for the interval 2T is shown in Figure lı. 

Assume that N electrons are emitted during the interval indicated. 


Then I(t) is the sum of the individual current pulses 


A 
£ 2(%-%-) 
Те = 7 b 


N=1,2,3 -.. 


and the time average 15 = 


~ Зо 1-2 


T го 
б! 

a. ЖЕ | T(t) At = ^ P Te) dt 

substituting 
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FIGURE 5 


Graph of i(t) versus Time 
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FIGURE 4 
Graph of i(t) versus Time Showing the Interval 2T 
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Now consider the electron stream to be made up of a great number of 


electrons whose fluctuations are considered to be ‚же: жш A count of 
electrons arriving at the plate in the interval 2T is N. Repeating this 
count for a total number of M times gives a value for the average number 


of electrons per second as 


лы Ne Ne + Ne + Nm T 
| Ут M — 00 MA (aT) 


As M is increased and T held constant, some of the N's will have 
the same value, thus the number of N's ofa particular value will in- 
crease as M is further increased. The probability of getting N electrons 


for some particular trial is defined as 


p(N) = lim Number of trials giving N electrons 1-5 
Mo A 4 o ons 


Since p(N) also depends on T it can be assumed that the randomness of 
the electron stream is such that the probability that an electron will 


arrive at the anode in the interval (t, t+ At) is 
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v At where At is such that vAt<«], and that this probability is in- 


dependent of what happened before time t or what will happen after time 
ttöt. For a very large number of trials this is a description of the 


Poisson distribution and is expressed as 


М. 
Thus a large number, M, of intervals of length 2T having exactly 
N arrivals will be to a first approximation Mp(N). For a fixed value of 


t, equation 1-3 may be written 


=,. 
= 


T, > И 
+: (t-t,)dty 1-7 


Nr, tt 


Miz 


yt 


and for A<t< 774 it effectively becomes 


IO = EN L(t) de 


It is now possible to write the average I(t) over all of the M in- 


tervals instead of only over those having N arrivals, so as M 


со 


It) = ‚С эу 1-9 
N= o 
оо N = 
3 5 ~ aT 6 12045 10 
М0 > „© 
00 : 
2 V I(t) dt 1-11 
-09 
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Campbell's theorem? states that the mean square value of the fluctuation 


about the average 5. Е 


A _ ор 
I) - Ton = v V it) a 1-12 
- ов 
= 1¢) - 2I(t) IE) - I(t) ПЭ 
| " N 
p (t) Y 7 i(t - tp) i(t - tn) 1-14 


Wei mm=\ 


Averaged over all values of ty, to, دوا‎ ecossos ty , With t held 


fixed 
Шз „кю Ыы T 4t pe | 
Lm 2) 2 | Tr ој == C(e- t) «(€ - $5.) 1-15 


n=l wei 


The multiple integral has two different values, so when ns m 


я А 22 d.t 
a Ts a Ё 
y ( m Zr 1-16 


and for n £ m its value is 


T 
| 265) dx 
т 2,7 





1-17 


T 
| ((е-4һ) 4%» 
= 2T 


Counting up the number of terms in the double sum shows that there are 


N of them having the first value and N? - N having the second value. 
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For - Tb L + < Т-А 


pepe = mat + NÜN-D “(044 И 
N "x E \ ( 1-18 


Averaging over all the intervals instead of only those having N 


arrivals, the final form can be written as 


Tia = x Poy Tye) 1-19 


N= O 
оо سے‎ 
=V\ ¢7(é)dé + I , 
1252 
- 09 


where the summation with respect to N is performed as in (1-11) and after 
the summation, the value (1-11) for I(t) is used. 

The relationship of the above development to the problem can be shown 
for a small increment of surface area located at one of the -points where 
the metallic airframe and plastic housing are adjacent. 

By equation 1-1 the individual current pulses can be described for a 
Single element, and when all of the points are taken into consideration, 
plus all of the locations where plastic Ном цав are adjacent to metallic 
airframe, equation 1-20 expresses the total average squared current. 


(Appendix I) 
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Ң. Previous Statistical Noise M¿asurements 

Some investigation of other statistical approaches to the measure- 
ment of noise was made before any design considerations were conceived. 

In their article, "A Statistical Approach to the Measurement of At- 
mospheric Noise" by Re Se Hoff and Re Co Johnson, +” it states in part,"... 
that atmospheric noise, causing noise voltage to be induced into a receiv- 
ing antenna...the envelope of this noise voltage after being amplified 
by the receiver stages and further modified by limiting and other non- 
linear processes, is an irregular function of time or time series Y(t), 
and may be treated by statistical processes." Also, den” has shown 
that the distribution of noise voltage when passed through a narrow band 
filter is normal (Gaussian),(Appendix I) that is, the fluctuation noise 
output of a band pass amplifier has a distribution of voltage versus time 
which follows the normal error law. The method used is similar to that used 
by Lord Rayleigh? in determining the probable distribution with time of 
the instantaneous amplitude of a large number of sine waves, all with 
equal amplitudes and random relative phases. From Knudzton's report“, "TH 
that the narrow band noise has the character of an amplitude and frequency 
modulated wave for the ideal rectangular band pass filter, where the 
average period of the envelope is a linear function of the bandwidth only." 

Thus it is possible to list the conditions that will effect the de- 
sign characteristics of the device to examine the properties of the pre- 
cipitation noise. 

(1) If a quantitative value of the noise power is desired, a narrow 
band pass filter of approximately 3 kc must be used since the center 
frequency is 126 kc. (See Appendix II) 

(2) The device must simulate a navigational system which has the fol- 

16 





lowing characteristics: 


(a) Bandwidth 30 ke 
(b) Resonant Frequency 126 Кс 
(с) Servo time constant response - 100 seconds 
(3) Ideally considered, the noise has a white spectrum. This of course 
is not true, but is permissible for tnis problem since Plank's radiation 
formula indicates that the thermal noise * is down l/e of its maxi- 
mum value of KT at about 10/ > 
To satisfactorily solve this dilemma of having two bandwidths it would 
be necessary to have two distinct systems; the first to heasure a quanti- 
tative value of the noise, and the second to ride next to the navigational 
system recording the relative noise input level for a specified naviga- 
tional track. It was decided to use the latter system to dy tne effects 
peak disturbance periods as a function of navigational error and to use 
the first system to set up a method to investigate the power spectrum of 
high density precipitation static by measuring the average power level 


across the band in small increments of frequency. 


17 











= 


Chapter II 


Description of Designed Measuring Devices 


l. SYSTEM I - MEAN SQUARE NOISE METER 
As shown in Appendix Il, it is possible to represent the output of 


a narrow band filter by 


e(t) = e (t) cos w,t + e 


UM (2111 Wat 


sin 

Using a very narrow band pass amplifier as the input stage, the mean 
value of the output noise voltage can be found by integrating the voltage 
over a long period of time. This integrator has an output which is made 
up of a component that increases igen with time, plus some variation 
due to the input a-c components wnich are smoothed out after sufficient 
elapsed time. It is specifically desired to — the mean square value 
of tne noise voltage ( е2), by constructing a device which will square the 
input and then integrate the low frequency (difference terms) components 
over a long period of time. The advantages are: 

(a) The recorded measurements are proportional to energy. 

(b) From the basic premise, that the noise has a normal distribution, 
it must then consist of a uniform spectrum of Spectral Density No 
volts? /cps at all frequencies other than zero, and have a true mean value 
(e)volts, it is necessary to continue the integration proce:s for a length 
of time 

Ww = (2 х 104 М.) /ғ% е ЁС 

in order to insure a 95% certainty that the measured value lies within 


P percent of the true nean. 
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Frequency Analysis of the Square law Circuit: The current of a 


diode may be represented by Ше infinite series 


: = 2 3 п 
% сет nu o horroroso cre 


in which the constants c, may be determined as are the coefficients ina 
Taylor series. 

Since the diodes are operated in a region where the dynamic character- 
istic is not a straight line, let it be assumed that the first two terms 
accurately represent the portion of the characteristic used. Thus 

i 2c, e+ c, ef 
p 4 2 
Substituting equation III-7: Appendix II 
is = c ( e.(t) cos wat * e&t) sin wt) 
+ €, ( e(t) cos w,t-ke.(t) sin Wat )2 

The only terms of interest are those containing difference terms since 
the integrator is essentially a low pass filter, therefore in this case 
iuris 

2 с2 е (+) eslt) sin w, t cos Wat 
By trigonometric identities tnis is equal to 

2 2 ес (+) e (t) E GA = t+cos W„t sin wot | 

Thus the presence of the difference terms allows the use of a large 
time constant in the integrator circuit to average the input low fre- 


quency terns. 
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Description of Circuits in System I 


Figure 1-0 Block Diagram of the System 

Figure 1-1 Narrow Band Pass Amplifier - Essentially a tuned grid, 
tuned plate tank, measured to be 3kc wide at half power points, with center 
frequency at 126 kc, and tunable over a 30 kc band. 

Figure I-2 - Logarithmic Amplifier - The amplifier is an instan- 
taneous linear amplifier designed to obtain a large linear range. A feed- 
back loop from the Miller integrator circuit which has a time constant of 
100 seconds causes a compressing effect to the input Signal. This allows 
examination of signals whose magnitude may vary from 10 micro-volts to 
104 micro-volts. The variable gain is calibrated to give approximately 
10 db gain for - 6 wolts of ave and 70 db gain for - 1 volt avc. 

Figure I-3 Square Law Circuit - This circuit was originally designed 
using the relation 


_ 2 
Ê, 2 КЕ; 


where Е, is the output voltage and Е2 15 the input voltage of the circuit. 
Knowing the desired avc characteristic and the time constant of the Miller 
integrator it is possible to approximate the value of k. This was accon- 
plished by plotting several curves using values close to the approximated 
value of k, constructing a circuit with variable resistors and adjustable 
bias levels, empirically determining the component values after the input 
and output values 'fit' the desired curve. 
Figure l-4 Detector and Miller Integrator - The detector is a peak 


reading device with a time constant determined by the center frequency 


T C detector = 1/1.26 x 102 = 8 micro-seconds 
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The basic principle of the Miller integrator in this system is the 


operation of averaging. It is performed by passing the general input through 


a low pass filter (integrator), thus the contributions to the output will 
be made only by those terms which have frequencies which fall within the 
range of the filter. In this problem the observation interval is fixed 
at 100 seconds, so thus the output, M(T), of the averaging filter may be 
related to its input by means of the convolution integral involving the 


filter's impulse response and weighting function, h(t), 


Т 
M(T) = h(t) Е (Т- +) а 


о 
where T is the duration of the observation interval. 
From "Statistical Errors in Но. Table I; the effective 


bandwidth, B, and the weighting function, h(t), are given for an integrator 


as 
h(t) = 0.0 where O0 > t ?T 
B, = 1/2T = 0.005 
The input will be the difference term 
E(t) = 2 Co e (t) e (t) sin (w, - w,)t 
for 


0<ч< 0.01 


Simplifying M(T) 


x 
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The above represents the output of the averaging filter for an 
observation interval of 100 seconds. 
Figure 1-5 Meter Coupling Circuit - This circuit uses an Esterline 


Angus type recorder, Model AW. The meter movement is a l ma type. 
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2. System II - NAVIGATIONAL DEVICE SIMULATOR & RECORDER 


As was previously stated, this system must duplicate a navigational 
device that is represented by the diagram in Figure 1. It must record 
noise voltage but must disregard the signal, so that if the circuit shown 
in Figure 1 were used in its entirety, a recording of the output would 
simply be the average value of the input signal and noise lumped together 
with no means of distinguisning one from the other. It is required then, 
that the signal must be eliminated and only the noise voltage recorded. 
Filtering of course will be useless since it removes those noise frequencies. 
that are of interest, so that the only solution is to phase out the signal. 
This will require two separate channels and the use of phase shifters, but 
in keeping with the basic premise, to duplicate the navigational system 
as much as possible, the only change of significance will be the intro- 
duction of a balanced modulator. lts purpose is to ict as a gate, opening 
and closing at a set rate, but 90° out of phase with the gate in the other 
channel. Thus when the two channels are added together, the signal will 
have been phased out but the noise components will be the sum of all the 
cosine terms (from one channel) and all of the sine terms (from the other 


channel). 
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Description of Circuits in System II 


. 
Figure 11-0 Block Diagram of System 


Figure 11-1 Band Pass Filter and Tuned Amplifier - Essentially tnis 
circuit is a single tuned, link coupled circuit, acting as a band pass 
amplifier with bandwidth of 30 kc with a resonant frequency of 126 kc. 
Figure II-2 Logarithmic Amplifier - Same as System I. 
Figure II-3 126 kc Crystal Oscillator and Phasing Networks. А highly 
stable Pierce crystal controlled oscillator is used as a local oscillator, 
the output of which is fed to two RC type phasing networks. If the out- 
put of the oscillator is 

еы = Е sin wt 
then the outputs of the phasing networks will be respectively 


e, = E sin w t =- 15° 
and 
ер 5 E sin wt t 15° 


Since ey and e, are ninety degrees apart and bear no fixed relationship to 


the input signal as yet, they can be represented as 


е1. Е ѕіп Wot 


eo 7 E cos Wot 

Each are amplified in separate channels and remain in these channels un- 
til detections take place. 

Figure II-l 126 kc Amplifier, Gate and Integrator - The first two stages 
are voltage gain stages which split ey into two voltages 180° apart in 
phase and thus drive a push pull cathode follower stage. By this means 


it is possible to offer a low impedance driving source to the balanced 


modulator (gate). 
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Frequency Analysis of Coherent Detector (Balanced Modulator) 
If the n noise voltage has the form at the output of the log 
amplifier of 
e(t) = Ay (t) sin wt + Ay (t) cos wt 
where 
ы = 21 (126 кс f 30 kc) 
the output of the balanced modulator can be computed in the following 
manner e 
Assuming matched diodes with similar current characteristics, the 


current can be represented by the infinite series 


1, = à) + ауе + а, е2 + a, 24 ....+а ей 
Since the integrator will act as a low pass filter on the difference 
terms of the output of the non-linear device will be of consequence. 
The d-c term (a,) can be ignored since transformer coupling is used. 
The voltage at point (1), Figure II-5 is 
E = eye e(t) 
The voltage at point (2), Figure II-5 is 
Eo = = e)+e(t) 


thus the value of Eat can be written as 


а р ® c 2 
о en = AE? = e) - Ae, - E) 
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Setting NR S TR 





. 
оо COE CEI 2 
O о тоог? - е(%)) 

Thus 

4 - ђ еј е(%) 

о ЕНЕ. 
K = h e 
Substituting 
е = -  EA,(t) sin wot sin wt - LEA,(t) sin w t cos wt 





K-LE sin wot 


Since only the difference terms are of consequence 
- біп м t cos wt = - sin СА - w) t - cos wot sin wt 
and 
sin wot sin wt = - cos (мо - Ww) t+ cos w, cos wt 
Thus eg can be represented as a voltage having frequencies made up of 
the difference terms. 


For channel two the analysis is the same with the exception of the 


value of е; 


е -ђ ё2 е(%) 


o 


к= Ц е, 
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Substituting 


e s “LE Ay (%) cos wot sin wt - hl E A,(t) cos wot cos wt 
О Ec da с ME ON Ao o oo ng re ħŮÁ 
К - ЦЕ соз wot 


Again only the difference terms are of consequence . 
- cos W,t sin wt s sin (м. - w)t - sin w t cos wt 


and 

- соз м cos wt = sin Wot sin wt - cos (w, - м) + 
as before, the voltage applied to the integrator can be represented as 
having бак баст» made up of the difference terms. 

The Miller integrator is similar to that 2 in System I except 
that the output can be either a positive or negative value, since the 
average will vary about zero and will have a mean average level of zero. 
Figure II-5 00 Cycle Oscillator and Amplifiers - The 400 cycle os- 
cillator is a Colpitts type with output fed to the two channels after 
being amplified and shifted in phase 90° with respect to each other by 
the two RC type networks. The method is similar to that previously men- 
tioned in connection with the 126 kc oscillator and the shift is in keeping 
with the previous phase shifts. The purse is to obtain the sum of the 
two channels across a transformer. It should be mentioned that the shift 
to 490 cycles is to enable a power gain to be obtained before taking the 
sum of the two channels. 

Figure II-6 400 Cycle Chopper and Amplification Stages. - The 00 cycle 


chopper is a Model manufactured by the Bristol Co. of Waterbury, Conn.. 
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юэ” ни 03 um. ы NSIS >» ga, ES = 5 
ihe necessary power amplification is provided by a single stare in each 


Channel with the two channels being coupled together using a push-pull 


amplifier. The output to the detector is the sum of the two channels. 


Figure II-7 Detector and Recording Circuit - The detector is a peak 


reading type with a time constant proportional to the ДОО cycle fre- 


quency. The Recorder is an Esterline Angus, Model AW, having а 1 ma 


movement. 








СУР REI 


OAV 












M3LA4IHS 250 
43040934 u 
ao, Sd900* 
Р “429 
3998 4.39 RIO 8 'ану 
.6%4+0 Sv tg 
Eum 
Ju3l3l 1d NV 
901 
139 8844081 | “өзіні е 
dann 'dNY| (642009 43111M 3179 53 92i 
eS v^ 0 St +g 
YIL IHS 
E 3SVHA ONINNL 
9 3921 оз эг! 3 Ч 





(Y 
m 








d3l3di ldNV 0OINNL 8 11945 1пам 
| 5 Ol 


ЭЕ. 
8 P | T - 


NOI1VH8lITVO —-. | 











2, овон Е 
Euros LE m 
| | 
ipe Чші2 47021 | ZTIE ITIGE- 














d3ldl IldWNV 


SH 914 





907 








SS 


WLX 
93921 


38 


AYOMLIN 1ЗІН6 35УНеа 
Ч01У1119050 71У1Х INI 


9 





AQG2 





2 TANN WHS 


9-0 | 
WEES 


| TINNVHO 


6-0 





р-Л. 914 


JMOLVHO3INI H311IN 9 31v9 'aaldl1dWNV 27 981 


= 0095 уа ланү 901 
«e с М 





AN MOYA 





слуз | 





и зедоно 01 


104100 
we 


1082 





TUETES 








H3ddOH2 e -j 





о 
> 
Eos um. 
O 2 
о 








Ss “IT “914 2 
| а | 
SH3lJdI1dWV E. ТЕЛ | 
| 
8 | 2 
| ES 2 
BOLYIIMSO 31949 00% ux ` 
| *. шээж er 
ED T. 
2 =. зм! | 
| E" Ho яс”. " | " 
ДЕ 5 о’ <, р 
ызааон2 | | La 
01 3 07 Эс 








” 9 — TT 914 


Y 


639У16 NOILVOISMGWY 
8 М3440Н0 312А2 00% 


C 


NNOO AunaH31VM 





00 101588 - = = = 
ка би Ж Т = БЕ 
3j19A9 OOv cm нін 
TT ? 012: яз ET 
u 3ddOHO | го Т : 
O! 13 иг: ямаа? ин“ 
1 спе = = OAL 
| 13 NNYH3 
Se | 800 did d JOHO 
|, | e 3042 E 
43083 c 
3042 | 4200» 
> d 
40193.130 
OL Se 








3NO 
чанкунзо 


Ч 3440НОЭ 





6S —SH (ЧУІМІ 





А 062 











| 


301 v3 


or A» 999 





10v zi 4 T дома 3 


49 SH QVIM1 
ло! »9! 
А 0682 











3. OPERATION 


The technique of taking data is quite simple. 

Previous calibration of the meter ged circuits should be made in 
the laboratory using a General Radio Signal Generator 8050, fixed at 12692 
with zero modulation. To simulate an antenna load a value of capacitance 
equal to the antenna capacitance of the particular aircraft model should 
be shunted across the input terminals so that the tuning capacitor can be 
set for the proper value. 

For an input of ten microvolts the meter should read a zero value 
while for a full scale reading a value of 100,000 microvolts can be used. 
Intermediate points can be located to suit the user. 

| Flights must be at least four hours long in order to accumulate suf- 
ficient data per flight. It has been found that a flight of this duration 
is necessary mainly to show the magnitude of the error of the navigational 
device for a known course and speed over a closed geographical loop of known 
dimensions. 

Flights in the vicinity of thunder clouds are desired with a record of 
the time spent in such locality marked on the recorder paper along side of 


the measurement record. 
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CHAPTER III 


RESULTS AND CONCLUSIONS 


l. Results: 

System Il was calibrated in the laboratory and will be delivered to 
the Air Force some time in April 1954. No particular recommendations or 
conclusions can be offered at this time since no flights have as yet been 
made, but as is typical with statistical data, the majority of time spent 
in investigation is used to gather data. 

System I was constructed and calibrated in the laboratory during the 
month of April 1954. A unique testing set up, developed by Tanner“? was used 
to generate corona pulses and is illustrated in Figure 5. 

Preliminary recordings of the output of this test setup indicated only 
minute amounts of power were being generated. It was decided to take the 
unit on a test flight to ascertain whether any general indications of pre- 
cipitation static could be obtained. 

A flight was made on April 14, 195 from 0830 to 1230 PST with route 
as follows: | 

From Monterey, California to Fresno, California, and then north-east 
to the Sierra Mountain chain, following this due north to Yosemite National 
Park. 

A long line of cumulus type clouds hung over the ridges offering good 
precipitation static conditions, so the system was placed in operation and 
immediately the recorder indicated moderate corona because all discharge 
wicks were removed from the plane previous to take-off. Figures III and 


IV are parts of the recordings taken during this flight and show vividly 


hh 
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the difference in the average power levels for an area of negligible pre- 


cipitation static and an area of heavy precipitation static. 


sica 

Investigations of precipitation static made using System II will 
record values related only to the average power input to a particular 
device so it is obvious that this information is of no value to anyone 
interested in the actual power spectrum or the possible frequency distri- 
butions of precipitation static. It was for this purpose that System I 
was designed. 

Since the data gathering period will take considerable time, only 
a brief outline of the method and approach will be given. 

The approach tothe problem can be again allied to the familiar treat- 
ment of fluctuation noise, using the autocorrelation function and its 
relation to the Fourier cosine integral. 

The autocorrelation function ø (T) of a time function i(t) is defined 
as the time average of the product of all pairs of points i(t) separated 
by an interval 


= ar MES C (t 4T) 4€ 


Substituting from Equation 1-20 
p) = VOS LG) c(€- T) dt f Te)” 


In the above the second term is the square of the d-c value I (t). 
If a current i(t) is defined as the fluctuation of the current about its 


mean value, then 
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Му -S 
8; (т) = Luo B (tr?) 
#00) = 12 = у (17 (0) а 


- 09 
Note that the total noise power depends not only upon the average current 


but on the shape of the current pulse as well. 

The noise power spectrum G; (w) can be obtained from the autocorre- 
lation function by the Fourier cosine integral relation, and is found to 
be 
G. (w) =2 y he (w т) 2(1 - cos wT,- w sin w Ta) 

( ч Ta)! 
where e is equal to the charge of electron, and can be plotted as shown 
in Figure 6. 
The noise power in the frequency range fa to f, is 
atn ағ 
5, 

In the frequency range for which wT.<<WY , G(f) is nearly constant 
in value. This is true no matter what the shape of the current pulse is 
Since the value of the Fourier Integral at zero frequency depends only 
upon the zero moment of the function. Thus C(0) depends only upon the 
area under 7, (*) and not upon its shape, so that the noise power іп а 


bandwidth Af jis 2veA fo 
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FIGURE 6 


Plot of Noise Power Spectrum versus WT 














A device then, having a bandwidth that is less than one radian can di- 
rectly measure increments of the power spectrum. This work is now in pro- 
gress using System I using the noise model described in the beginning of 
this chapter. It is obvious that the noise model device is inadequate as 


a generator for heavy corona conditions and investigation and design of a 


superior device is contemplated. 
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APPENDIX I 





THE NORMAL DISTRIBUTION 
l. Probability: 
(a) If an event is certain to occur it is assigned an ar- 
bitrary value of unity. For example: 
In the toss of a single die the probability of 


a six ів 1/6 


a five " 1/6 
a four " 1/6 
a three " 1/6 
a two n 1/6 
a one " 1/6 


thus, if six probabilities are mutually exclusive, the probability that 
either one or the others will occur is equal to the sum of the separate 
probabilities. 

(b) When a particular result in any trial is independent of 
the results obtained in any other trials, is called a random process. For 
instance: 

A large number of independent trials in tossing a coin might reveal a 
consecutive string of heads. This will not affect the next toss since the 
probability is still 1/2 that the coin will cme up either a head or a tail. 

(c) The probability of the simultaneous occurrence of n events 
is the product of the n individual probabilities of n events. Thus the 
probability of winning all tosses in a season of eight footballs games 
would be (0.5)° = 0.0039 

The above three rules form the basis of the probability distribution 


laws upon which the mathematical description of noise is based. 
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2. Distribution formulas: 

If the symbol ГЕ is made to stand far the probability that an 
event will occur exactly n times in m independent trials, then the basic 
distribution formula, (Bernoulli Distribution), is 


=- M n m-n n m-n 
PR 5 0 р(1-р) р (1-р) 


" 
8 


where p is the probability for the event to occur in any particular trial. 


This specializes to, (Poisson Distribution), 


8 - - for m >> 
Pym) = a (Mm? EO 
п] и >> 1 
where 
n = mp 


When n becomes very large as in the case of high density Poisson 


. noise, the above further specializes to the Gaussian Distribution. 


22 = 
е) = 1 y) an 
when n >? 1 
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3. Тһе Central Limit Theorem 


The Central Limit Theorem states that subject to very lenient 
conditions, and regardless of the amplitude distributions of the indi- 
vidual components, as the number of components n approaches infinity, 
the amplitude distribution of the current approaches the Gaussian or 
normal distribution. 

Thus it is possible to relate the limiting form of the Gaussian 
Distribution to the Mathematical Description of Noise (ChapterI) by 
virtue of the large number of independent events (electron crossings) 
which add up to give the independent MT I. This distribution 
will be independent of the amplitude distributions of the individual 
events (shapes of the electron pulses). It is required, generally speak- 
ing, that the contributions of the individual components be small compared 


with the total value of the sum 
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è APPENDIX II 


PROPERTIES OF NARROW BAND NOISE 


1. Examining the statistical properties of the output of a narrow 
band filter plotted in Figure A-I, and the assumption, that the band 


width of the filter, f,<“ 5. „the output сап be expressed in the form 


e(t) = E(t) cos (wt t g(t)) 11-1 


provided that the input noise has a Gaussian distribution. 
E(t) then represents the magnitude of the envelope of the output and 
g(t) the phases changes. 


Using the Fourier representation of the output as 


e(t)=> (a cos пи Ф + bı sin nw t) 


where the summation is to include all frequencies which is in the range 
of the filter. 
Setting up the sum and difference terms so as to show relationship 


to a sine wave of frequency Wa? 


cos nw t 2 соз ( (ач, -Wat te не \ 11-3 


cos (nw. -Wp)t cos wit - sin (mw, w) sin wt 


sin nw t sin | сн, ut + 23 II-h 


sin (nw, - wit cos wit + cos (nw, - 4,0% sin wt 








D 
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Frequency 


FIGURE A-I 


Output of a Narrow Band Filter 








Letting ОР" (t) represent the summation of all the coefficients of 


coS W,t in the expression for e(t), and e,,,(t) represent the summation 


of the coefficients of sin Wņpt» then 


A 
= E із, cos (nw, - w„)t+b, sin (nw, - “| ІІ-5 
^ 
and 
= 
int) = ? Us cos (nw, - w,)t - a, sin (nw, -wt | 11-6 
bat 
Substituting 


e(t) = Ti (t) cos wi ¥ e (%) sin wt 11-7 
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Thus, ё 05 (%) апа © in (t), are slowly varying compared with 


cos Wn t or sin wt. 


From the Expression III-1, the magnitude 


E(t) = 1 94 (t) + е (+) NICE 


(t) * tan Е v (t) / Баг СЕ) ІІ-9 








апа 


The random variables e (t) and e. 


сов in (t) are sums of independent 


i 
variables with Gaussian amplitude distributions. The average power can 


then be written as 
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